Electronic coherence in (5-Pu: A DMFT study 
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A combination of Density Functional Theory and the Dynamical Mean Field theory (DMFT) 
is used to calculate the magnetic susceptibility, heat capacity, and the temperature dependence of 
the valence band photoemission spectra. The continuous-time hybridization expansion quantum 
Monte-Carlo is utilized to provide the first approximation-free DMFT solution of fee 5-Pu which 
includes the full rotationally-invariant exchange interaction. We predict that S-Pu has a Pauli- 
like magnetic susceptibility near ambient temperature, as in experiment, indicating that electronic 
coherence causes the absence of local moments. Additionally, we show that volume expansion causes 
a crossover from incoherent to coherent electronic behavior at increasingly lower temperatures. 
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Pu embodies a forefront of both technology and theo- 
retical condensed-matter physics. Elemental Pu displays 
exotic physical behavior that continues to defy explana- 
tion. For example, it exhibits six allotropic phases at 
ambient pressure, the low-density fee S phase has a neg- 
ative coefficient of thermal expansion, and the volume 
expands by more than 25% when the system is heated 
from the high-density monoclinic a phase to the 6 phase. 
Even many of the basic electronic properties have not 
been definitively characterized and understood. Regard- 
ing the pure d phase, a complicating factor is that it is 
only stable in the temperature range 580K < T < TOOK. 
However, the 5 phase can be stabilized at low temper- 
atures by a variety of alloying elements such as Ga and 
Am. This allows for the experimental exploration of the 6 
phase at low temperatures, with the caveat that it is not 
clear what changes the alloying element may be inducing. 

Lashlcy et al. [l[ have measured the magnetic sus- 
ceptibility to be Pauli-likc in both the a and S phases, 
and hence detect no presence of localized magnetic mo- 
ments. Similarly, Heffner et al. 0, i, 13 have used fiSR 
and showed that there is no ordered magnetic moments 
in a Pu nor in (5-stabilized Pu (i.e. 4.3% Ga) for temper- 
atures down to 4K. Nuclear magnetic resonance (NMR) 
measurements by Curro et al. also show an absence of 
magnetic moments. 

The linear coefficient of the specific heat for S- 
stabilizcd Pu has been measured by various groups and 
the resulting values are 

42^74? for alloying with 2% of 



GaQ, 64- 



m J 



for 5% of AM, and 35 - 55- 



for 



8-20% of Am[7[. The large variation among these mea- 
surements may be due to the fact that the 5 phase has 
been stabilized by a different alloying element in each 
study. Furthermore, Pu easily incorporates elements 
such as Fe during synthesis, while radiation self-damage 
may generate further defects that may influence the low- 
temperature specific heat. 

Several previous studies have applied a combination 
of Density Functional Theory and the Dynamical Mean 
Field Theory (DFT+DMFT) [| to (5 Pu. DMFT re- 



quires a solution of an auxiliary quantum impurity prob- 
lem, and for the corresponding impurity model of Pu, no 
exact method was available in the past. Savrasov et al.^ 
used an interpolative solver to calculate the energy and 
the photoemission spectra of Pu. The approach yielded 
a significant improvement for the volume of the 5 phase 
of Pu compared to DFT. Pourovskii et al. com- 
puted the photoemission spectra and the heat capacity 
using the Fluctuation Exchange Approximation (FLEX) 
as an impurity solver. The applicability of FLEX to Pu 
is questionable given the strongly correlated nature of 5- 



Pu. Zhu et aL[ll| computed the photoemission spectra 



using the the Hirsch-Fye quantum Monte-Carlo impurity 
solver and show that the occupation of the / orbital is 
close to n/ ~ 5. However, the limitation of this method 
to treat the realistic atomic multiplet structure (Hund's 
rule coupling) Q precludes a precise description of the 
problem. Shim et al. predicted the photoemission 
spectra, the X-ray absorption spectroscopy branching ra- 
tio, and the mixed-valence nature of Pu. Although the 



mixed- valence state was identified in ref. |12| . the temper- 
ature and pressure dependence of the electronic state was 
not addressed. 

In this letter, we demonstrate the absence of magnetic 
moments in 5-Vn by computing the magnetic susceptibil- 
ity as a function of temperature. We show that expand- 
ing the Pu lattice results in an incoherent metallic state 
with Curie- Weiss susceptibility at increasingly lower tem- 
peratures. Additionally, we elucidate the nature of the 
the mixed-valence state by predicting the temperature 
dependence of the photoemission spectra in an extended 
temperature regime starting below ambient and contin- 
uing beyond the melting temperature. 

DMFT maps the interacting lattice problem onto an 
impurity problem where the non-interacting bath func- 
tion is determined self-consistently The effective 
impurity problem is then solved using the continuous- 
time quantum Monte-Carlo (CTQMC) method [3, [H]. 
More specifically, the recently developed hybridization 
expansion CTQMC method is used to exactly sum the 
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diagrams resulting from expansion in powers of the hy- 
bridization streiigth between the Pu atom and the DMFT 
fermionic bath |14l |. This method aUows one to include 
the full rotationally-invariant exchange interaction with- 
out approximation, in contrast to Hirsch-Fye QMCQ. 
Although our implementation of CTQMC is extremely 
efficient, massive parallel computer resources are required 
to solve the 14 orbital DMFT impurity problem for Pu. 
The Atlas supercomputer at Lawrence Livermore Na- 
tional Lab was used to perform the calculations, and time 
was awarded under the Atlas grand challenge program. 

DFT+DMFT calculations were performed using an or- 
thogonalized LMTO basis which provides / orbitals that 
have a maximal / character [l^. In accordance with pre- 
vious studies, the on-site Coulomb repulsion was chosen 
to be J7 = AeV and, the Slater integrals {F2, F4, and Fg) 



were computed using an atomic physics code jl7l | and 
rescalcd by 80% to account for screening in the solid. 
Summations over the first Brillioun zone were performed 
with 15x15x15 meshes, with the exception of the heat 
capacity calculations where up to 40x40x40 meshes were 
used to ensure convergence. Pu valencies of Nf = 4 — 7 
were retained in the QMC simulation, and retaining 
higher/lower valencies had no appreciable influence on 
the results. The total occupation of the Pu-5/ orbitals 
in this study is approximately nj ~ 5.2, consistent with 
previous calculations [13 and electron en ergy loss exper- 
iments of the N edge[l8| and the edgejigj. 

We proceed by first exploring the qualitative effect of 
electronic correlations on the local spectra, both with 
and without the rotationally-invariant exchange interac- 
tion. The real-frequency spectral function is obtained 
by using the maximum entropy method to analytically 
continue the imaginary time Green's function measured 
in the CTQMC simulation. The LDA spectrum dis- 
plays a strong spin-orbit splitting among the S ~ ^ and 
S" = I states (see figure [TJi). The DFT-I-DMFT spec- 
trum without exchange indicates that spectral weight 
from low energies (ie. near the Fermi energy) has trans- 
ferred to higher energies (see figure [T]d) . Including the 
full exchange interaction reduces the spectral weight in 
the 5* = I quasiparticle peak and hence the electronic 
coherence scale. Additionally, the 5=1 peak above the 
Fermi energy is broadened due to a multiplet splitting. 
The spectrum with full exchange interaction is similar to 
the three-peaked spectrum obtained by Shim et aZ.|12|. 
The = f has a central peak just below the Fermi en- 
ergy and a peak near — leV, while the S — ^ states have 
a peak near — 0.6eV. These features are much broader 
than those obtained by Shim et a/.fl^. but this is ex- 
pected due to the limitations of the maximum entropy 
method in obtaining the real-frequency data. 

The self-energy Y,{iuj) offers further insight into the 
nature of the electronic correlations in (5-Pu. The quasi- 
particle weight is determined by the slope of the imag- 
inary part of the self-energy (ie. Z = 1/(1 — ^^7^)), 
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FIG. 1: The spectral functions for the Pu /-electrons within 
DFT (top), DFT-I-DMFT without exchange (middle), and 
DMFT-I-DFT with the full rotationally-invariant exchange 
(bottom). 



and within DMFT the effective mass of the electrons due 
to electronic correlations is the inverse of the quasipar- 
ticle weight (ie. ^ — 'k)- ^^^^ absence of exchange, 
the 5* = I states are very weakly correlated, having an 
average Z « 0.7 while the -S* = | states have an aver- 
age Z w 0.41 for the equilibrium volume of S Pu (see 
figured^). Weak correlations for the S = 7/2 states is 
expected given that in the absence of Hund's coupling 



the j-j coupling scheme is adequate, and hence the S = 
states are nearly empty. Alternatively, the S = ^ states 
are much closer to a non-zero integer filling and are mod- 
erately correlated. When the exchange is included, the 
quasiparticle weight is substantially decreased, resulting 
in Z^/2 = 0.26, Zj^2 ~ 0.32 for the equilibrium volume. 
The exchange interaction pushes Pu towards intermedi- 
ate coupling, where the S = 7/2 states are more mixed 
into the ground state and hence they becomes heavier. 
With increasing volume, the kinetic energy decreases and 
hence correlations increase resulting in a smaller Z . The 
only exception is the 7/2 orbital in the absence of ex- 
change, where the occupation slightly decreases with in- 
creasing volume. 

Having established the effect of correlations on the 
spectra and the importance of the full rotationally- 
invariant exchange interaction, we proceed to explore the 
heat capacity and the magnetic susceptibility. Within 
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FIG. 2; The quasiparticle weight Z (top), and hnear coeffi- 
cient of the heat capacity (bottom) as a function of volume. 
The heat capacity was calculated with the full exchange in- 
teraction. 



Fermi liquid theory, the linear coefficient of the heat ca- 
pacity is given by 7 27rfc^ p^(S) ^ ^^^^^ ^ runs over 
all orbitals, pa is the local density of states, and is 
the corresponding quasiparticle weight. The heat capac- 
ity as a function of volume is shown in Fig. [21 As the 
volume is increased, the heat capacity increases due to 
the fact that the quasiparticle renormalization amplitude 
Z decreases and the spectral density at the Fermi energy 
increases. The predicted heat capacity for the equilib- 
rium volume of (5-Pu is 20.4 The difference bc- 



tween the predicted value of the heat capacity and the 
experimentally measured values of 35 — 55 y^o"'_5j^2 may 
be be due to several factors. Given that the /-electron 
spectral function is extremely steep in the vicinity of the 
Fermi level (see Figure [1]), it is clear that the value of 7 
is sensitive to small changes in the Fermi energy. There- 
fore, approximations in the calculation, such as the lo- 
cal density approximation to the Kohn-Sham potential, 
may have a non-negligible influence on the heat capac- 
ity. The second potential cause of this difference might 
be that the electron-phonon coupling may further renor- 
malize the hoppings and this is not included in our cal- 
culation. Given that our prediction is smaller than ex- 
periment, inclusion of electron-phonon coupling would 
improve agreement with current experiments. It should 
be noted that the largest volume point is not yet in the 
Fermi liquid regime, as demonstrated below, and hence 
the Fermi liquid formula for the heat capacity should only 
be considered as an estimate for the largest volume. 
The local magnetic susceptibility is calculated as a 
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FIG. 3: Local magnetic susceptibility as a function of tem- 
perature for different volumes. 



function of temperature for different volumes using the 
following expression: 

Xlocal = Jo^dT{M{T)M{0)) 

where M = L + 2S, L is the total orbital angular mo- 
mentum, and S is the total spin angular momentum. For 
the equilibrium volume of Pu, the susceptibility is rela- 
tively flat (ie. Pauli-like) below 60QK and diminishes at 
higher temperatures. This behavior is consistent with ex- 
perimental measurements. Our calculations predict that 
the reason that magnetic moments are not seen in i5 Pu 
is because the system is coherent and the moments are 
therefore screened. This predicted behavior is in agree- 
ment with experimental observations P, [3, As the 
volume is increased by 16%, the low-temperature Pauli- 
like contribution has been renormalized to higher values 
and transitions to a Curie-like behavior around 400K. 
The increase in the Pauli contribution is consistent with 
the decrease in Z as the volume is increased, and this 
illustrates an enhancement of the electronic correlations. 
As the volume is further increased by an additional 16%, 
the susceptibility is Curie-like to the lowest temperatures 
reached in this study. These results clearly indicate a de- 
coherence of the electrons. As the volume is increased, 
the coherence energy of the electrons decreases and there- 
fore a transition from coherent to incoherent behavior (ie. 
Pauli-like to Curie-like) occurs at increasingly smaller 
temperatures. This signifies that a local magnetic mo- 
ment has emerged at progressively lower temperatures 
as the volume is increased. 

The general behavior observed in these magnetic sus- 
ceptibility calculations is consistent with experiments. 
First, doping Pu with Americium causes the Pu lattice 
to expand, effectively increasing the volume. The Pauli 
contribution of the Pu atom is shown to increase as the 
Americium content is increased and the effective volume 
is increased |20f. This is consistent with the enhancement 
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of the Pauli term that wc observed. Secondly, when hy- 
drogen is doped into the system, the fee Pu lattice ex- 
pands by more than 50% and Curie-like behavior is mea- 
sured for the magnetic susceptibility down to tempera- 
tures of 50K 2l| . This is qualitatively consistent with 
the largest volume expansion in our calculations. 




FIG. 4: The temperature dependence of the S = | /-electron 
spectral function using the rotationally-invariant exchange in- 
teraction. The inset displays the peak height as a function of 
temperature. The red points correspond to our data, while 
the blue curve is the parametrization from Ref. [^^ assuming 
an onset temperature of 800K. 



When the electronic behavior departs from the Fermi 
liquid theory and enters the incoherent regime, there is 
a clear signature in the spectra. In Fig. 31 we show the 
temperature dependence of the local spectral function at 
ambient pressure. At high temperatures, the spectrum 
is diffuse. As the temperature is decreased, a quasipar- 
ticle peak continually builds and eventually saturates at 
T ~ 500/'Sr. The coherence temperature may be defined 
as the temperature at which the quasiparticle peak nears 
saturation, and we define 75% saturation to be the on- 
set of coherence. Hence our estimation for the coher- 
ence temperature is approximately 800 K, consistent with 
ReflH. 

The inset of Fig. 2] shows the temperature depen- 
dence of the height of the quasiparticle peak. Notice 
that this behavior is very different from the temperature 
dependence in heavy-Fermion compounds (23| recently 
parametrized in Ref. [2^ . For comparison, we plot the 
best fit of the parametrization of Ref. [l^ to our data. 
The inability of this parametrization to describe our data 
is due to the mixed-valence nature of Plutonium. 

In conclusion, we have performed approximation-free 
DMFT calculations including the rotationally-invariant 
exchange interaction for S-Pu. The efficient CTQMC 
algorithm has allowed us to reach both high tempera- 
tures and temperatures below ambient. The quasiparti- 

states of (5-Pu is found to be 



Z ~ 0.25, indicating the presence of appreciable elec- 
tronic correlations. Calculation of the magnetic suscep- 
tibility indicates Pauli-like behavior for the equilibrium 
volume of (5-Pu, in support of experimental measure- 
ments. This indicates that in (5-Pu the moments are 
screened. Expanding the volume causes the electrons 
to crossover from coherent to incoherent behavior at in- 
creasingly lower temperatures. This crossover is illus- 
trated in the temperature dependence of the spectra. 
The prediction of incoherent electronic behavior, man- 
ifested as Curie-like behavior in the magnetic suscepti- 
bility, is consistent with experimental measurements in 
Pull2. The importance of the inclusion of the rotation- 
ally invariant exchange is illustrated in the spectrum and 
in the reduction of the quasiparticle renormalization am- 
plitude from Z = 0.41 to Z = 0.25. 
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